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Abstract 

We report an X-ray study of the evolved Galactic supernova remnant (SNR) G156. 2+5.7 based on six 
pointing observations with Suzaku. The remnant’s large extent (100' in diameter) allows us to investigate 
its radial structure in the northwestern and eastern directions from the apparent center. The X-ray 
spectra were well fit with a two-component 11011 -equilibrium ionization model representing the swept-up 
interstellar medium (ISM) and the metal-rich ejecta. We found prominent central concentrations of Si, S 
and Fe from the ejecta component; the lighter elements of O, Ne and Mg were distributed more uniformly. 

The temperature of the ISM component suggests a slow shock (610-960 km s” 1 ), hence the remnant’s age is 
estimated to be 7,000-15,000 yr, assuming its distance to be ~1.1 kpc. G156.2+5.7 has also been thought to 
emit hard, non-thermal X-rays, despite being considerably older than any other such remnant. In response 
to a recent discovery of a background cluster of galaxies (2XMM J045637.2+522411), we carefully excluded 
its contribution, and reexamined the origin of the hard X-ray emission. We found that the residual hard 
X-ray emission is consistent v/ith the expected level of the cosmic X-ray background. Thus, no robust 
evidence for the non-thermal emission was obtained from G 156. 2+5. 7. These results are consistent with 
the picture of an evolved SNR. 

Key words: ISM: abundances ISM: individual (G156. 2+5.7) — supernova remnants X-rays: 

ISM 


1. Introduction 

G156.2+5.7 (RX04591+5147) was initially discovered 
in X-ray with ROSAT (Pfeffermann et al. 1991). Its age 
was estimated at ~26000 yr, based on the Sedov model, 
hence is thought to be an evolved supernova remnant 
(SNR). The distance is comparatively nearby (~1.1 kpc; 
Pfeffermann et al. 1991) and its large apparent size (~100' 
in diameter) enables us to study the detailed plasma struc- 
ture. While G156.2+5.7 is one of the brightest SNRs in 
X-rays, its radio emission is quite faint; this is its most 
remarkable feature. Shortly after its discovery in X-rays, 
the radio continuum emission was also detected with the 
Effelsberg 100-m telescope (Reich et al. 1992). As they 
point out, the surface brightness at 1 GHz (5.8 xlO -23 
W m“ 2 Hz -1 sr -1 ) is the lowest among all the known 
SNRs. The radio morphology along the northwest and 
southeast rim shows limb brightening (Reich et al. 1992; 
Xu et al. 2007) that is typically seen in so-called ’’barrel- 
shaped” SNRs (Kesteven &; Caswell 1987). In contrast 
to its radio morphology, G156.2+5.7 has centrally-filled 
thermal emission in X-rays. While Lazendic et al. (2006) 


lists G 156.2+5. 7 as a mixed-morphology (MM) SNR for 
this reason, it does not fit the definition of MM SNRs in 
several ways: the X-ray shell is easily distinguishable and 
there are neither CO clouds nor OH (1720 MHz) masers 
in its vicinity. A recent optical observation suggests that 
the distance to G 156. 2+ 5. 7 is smaller than that previously 
estimated (Gerardy & Fesen 2007). Their main argument 
is that several bright Ha filaments correlate with nearby 
interstellar clouds (~0.3 kpc). However, the X-ray and 
the radio observations both agree that G 156.2+5.7 is lo- 
cated at 1-3 kpc, i.e., an evolved SNR (Pfeffermann et al. 
1991; Reich et al. 1992). The distance and the age of this 
SNR remain open questions. 

The X-ray spectrum of G 156.2+5. 7 is also largely typ- 
ical of an evolved SNR. Yamauchi et al. (1993) observed 
G156.2+5.7 with the Ginga LAC (large area proportional 
counters) and detected extended X-ray emission in the 1.2- 
3.0 keV energy band. Yamauchi et al. (1999) performed 
an additional observation with ASCA and found that its 
spectrum consists of two components with different fea- 
tures. The soft component has several emission lines, 
suggesting a thin thermal origin whereas the hard compo- 


2 


Author(s) in page-head 


[Vol. 



Fig. 1. Left : ROSAT PSPC image of G156.2+5.7 with the Galactic coordinate overlaid. The FoV of Suzaku. XIS are shown with 
solid squares. Middle : Radio (92 cm) image of G156.2+5.7 obtained by the Westerbork Northern Sky Survey (WENSS) with the 
Westerbcrk Synthesis Radio Telescope (WSRT). The coordinate scales axe matched in both images. Right: True-color XIS image of 
G156.2+5.7. The energy bands of red, green and blue correspond to 0.5-0.7 keV, 0.7-1.0 keV and 1.0-5.0 keV, respectively. The red 
circular regions were excluded from our spectral analyses. 


Tkble 1. Summary of the 6 Suzaku observations 


Name 

Obs. ID 

Obs. Date 

R.A., Dec. (J2000) 

Effective Exposure 

NW_rim 

501075010 

2007-Feb-16 

04 h 56 m 54 s .9, 52°24'12". 6 

50.5 ks 

Center 

501106010 

2007-Feb-17 

04 h 58 m 54 s .9, 51°44'36".l 

51.2 ks 

E_rim 

501074010 

2007-Feb-18 

05 h 03 m 13 s .5, 51°37'06".2 

53.3 ks 

East 

504082010 

2010-Feb-21 

05 h 01 m 03 s .0, 51°41'41".6 

50.3 ks 

NW 

504081010 

2010- Mar-03 

04 h 57 m 42 s .l, 52°04'31".l 

52.9 ks 

S_rim 

504080010 

2010- Mar-04 

04 h 59 ra 41 6 .8, 51°00'36".0 

52.6 ks 


nent is featureless. However, due to the lack of statistics, 
they were not able to decide whether the featureless hard 
component has a thermal origin or a non-thermal origin. 
Katsuda et al. (2009) observed two rims (northwest and 
east) and the center with Suzaku. They found two ther- 
mal components, one of which they ascribed to interstel- 
lar material heated by the forward shock, and the other 
reverse-shock heated ejecta. They also detected hard X- 
ray observed from the northwestern rim and the center, 
which they concluded is likely to be non-thermal syn- 
chrotron emission. If this is the case, G156.2+5.7 would 
be the oldest SNR producing non-thermal X-ray emission. 
However, a recent Suzaku observation discovered a cluster 
of galaxies (2XMM J045637.2+522411; Yamauchi et al. 
2010) located in the northwestern rim. Since the surface 
brightness of this extended source is higher than that of 
the rim, we need careful consideration of its contribution 
to determine the true origin of the hard component. 

Katsuda et al. (2009) also estimated the progenitor 
mass to be less than 15 M 0 from the metal abundances 
of the thermal emission detected in the central por- 
tion. While no compact source has been discovered in 
G156.2-I-5.7 (see Kaplan et al. 2006), their result suggests 
that its origin is a core-collapse supernova (SN) rather 
than a type la SN. More recently, Hudaverdi et al. (2010) 
observed the center of this remnant with XMM-Newton 
and measured its metallicity. From the results, thev sug- 


gested that G 156. 2+5. 7 is an O-rich SNR. 

We performed 6 pointing observations of G 156.2+5.7 
with Suzaku. The fields of view (FoV) cover three sec- 
tions of the rim and three inner regions, and include the 
ones discussed in Katsuda et al. (2009). In this paper, 
we report the ejecta distribution and review the origin of 
the hard X-ray emission detected by Katsuda et al. (2009) 
along with the newly observed regions. We also discuss 
the age of this SNR based on these analyses. 

2. Observations and Data Reduction 

We summarize the 6 observations in table 1. Data were 
taken in 2007 and 2010, during the Suzaku AOl and A04 
phases, respectively. The observed regions are shown in 
figure 1. The squares represent the FoV of the Suzaku 
XIS (X-ray Imaging Spectrometer; Koyama et al. 2007). 
The data were analyzed with version 6.9 of the HEAsoft 
tools. For data reduction we used version 9 of the Suzaku 
Software. The calibration database (CALDB) used was 
the one updated in March 2010. We used revision 2.2 
of the cleaned event data and combined the 3x3 and 5x5 
event files. In order to exclude background flare events, we 
obtained the good time intervals (GTIs) by including only 
times at which the count rates in the \arious detectors are 
less than +3*7 of the mean count rates. 
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Tkble 2. Surface brightnesses of hard emission in G 156.2-1-5.6 and two nearby observations 


Name 

l 

[deg] 

b 

[deg] 

surface brightness at 2.0-10.0 keV 
[10“ n erg cm^'Meg -2 ] 

G 156.2+5.7 NW-rim 

155.5 

5.8 

1.88+0.26 

G156.2+5.7 E-rim 

156.7 

6.1 

1.61+0.23 

G156.2+5.7 S-rim 

156.8 

5.3 

1.56+0.22 

IRAS 05262+4432 BGD 

165.0 

5.6 

1.78+0.25 

SGR 0501+4516 BGD 

161.4 

1.9 

1.89+0.26 



Fig. 2. Left : Surface brightness (SB) distribution of the CXB component. Filled small circle, open circle, cross, and square 
correspond to the results at E-rim, East, Center, and NW, respectively. Green (left) corresponds to the results at S_rim. Black and 
Red triangles (right) correspond to the results in the NW_rim, except in the regions of 1'.2 and 4'.0 around 2XMM J045637.2-j-522411, 
respectively. The '-ertical error bars correspond to 1(7 errors. Right : Histogram of the surface brightnesses. The horizontal arrow 
represents the 1 a spatial fluctuation of the CXB surface brightness estimated by Tawa (2008). 


3, Spectral Analysis 

For the following analysis, the energy range of 0.4-8.0 
keV was used for XIS1 (back-illuminated CCD; BI CCD) 
and XISO and XIS3 (front-illuminated CCDs; FI CCDs), 
since there was no signal above 8.0 keV after the back- 
ground subtraction (described below). In order to gener- 

ate response matrix files (RMFs) and ancillary response 
files (ARFs), we employed xisrmfgen and xissimarfgen 
(Ishisaki et al. 2007). We used XSPEC version 12.6.0 
(Arnaud 1996) for all the following spectral analyses. 
We excluded 2XMM J045637.2+522411 and three other 

point-like sources identified from the XIS image (0. 5-5.0 
keV). The extracted regions are shown in figure 1. 

3.1 . Background Subtraction 

Since G 156.2+5.7 is a large diffuse source and the FoV 
is filled with the SNR’s emission, it is hard to obtain 
the background spectra from inside the FoV. We also 
have no background data from the neighborhood of the 
G156.2+5.7. For the following analysis, we used a gen- 

erated non X-ray background (NXB) spectra by employ- 
ing xisnxbgen (Tawa et al. 2008). We also took into ac- 

count the contributions of the local hot bubble (LHB) 
and the cosmic X-ray background (CXB) by assuming 
literature-based models. For the LHB, we used an APEC 
model (Smith et al. 2001) with solar abundances (Anders 


&; Grevesse 1989) and an electron temperature of 0.1 
keV (Miller et al. 2008). For the CXB we employed a 
broken power-law model with photon indices of T=2.0 
(<0.7 keV) and T=1.4 (>0.7 keV) as shown by Miller 
et al. (2008). Hereafter, we call this model the CXB 
model. We also reviewed the effect of the Galactic Ridge 
X-ray emission (GRXE). The GRXE surface brightness 
at l = 100°, |6| < 1° is ~ 1.4 x 10“ 12 erg cm“ 2 s _1 deg” 2 
in the 0.7-10.0 keV band (Sugizaki et al. 2001). Since 
G156.2+5.7 is located outside the FoV of Sugizaki et al. 
(2001) and away from the galactic plane, the contribution 
of the GRXE is expected to be negligible. For exam- 
ple, GRXE surface brightness above 2.0 keV decreases 
about one-fortieth from 5 = 0° to 6^9° (Revnivtsev 
et al. 2006). Additionally, the average 0.7-10.0 keV 
band surface brightness of G156. 2+5.7 is ~ 4.3 x 10" 3 * * * * * 9 erg 
cm“ 2 s _1 deg~ 2 (Yamauchi et al. 1999) which is an order 
of magnitude higher than the minimum Suzaku surface 
brightness of ~3.4 x 10~ 10 erg cm _2 s _1 deg“ 2 that we see 
in the S_rim observation. In both cases, we conclude that 
the effect of the GRXE is negligible in G156.2+5.7. 

Table 2 shows the surface brightness of the CXB and the 
GRXE near G156.2+5.7 and within a few degrees using 
observations of IRAS 05262+4432 (narrow-line Seyfert 1; 
Kawaguchi et al. 2008) and SGR 0501+4516 (soft gamma 
repeater; Enoto et al. 2009). Since both are point-like 
sources, we excluded the source regions with radii 4'.0 
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and fitted :he spectra obtained from the remaining regions 
with the CXB model. Comparing the surface brightness 
of each observation, the level of the CXB and the GRXE 
is uniform within the range of error. Figure 2 left shows 
the surface brightness distribution of the CXB component 
inferred from the fit to each annular region in G156. 2+5.7 
(the details will be argued hereinafter). The filled cir- 
cles, open circles, crosses, and squares correspond to re- 
gions from the Ej-im, East, Center, and NW, respectively. 
Green represents regions from the S_rim. The result sug- 
gests that the level of the CXB (and the GRXE) inside 
this remnant is also unchanged within the range of error. 

3.2. Model Fit for the Rim Regions' Spectra 

Since the rims are mostly dominated by an interstellar 
medium (ISM) component (Katsuda et al. 2009), we ap- 
plied a single-component VNEI (NE1 ver.2.0; Borkowski 
et al. 2001) model with TBabs (Tiibingen-Boulder ISM 
absorption model; Wilms et al. 2000) in XSPEC for the 
rim spectra. In this model, the abundances relative to 
the solar values (Anders & Grevesse 1989) of N, 0, 
Ne, Mg, Si, S and Fe were allowed to vary freely and 
that of Ni were linked to Fe. The other element abun- 
dances were fixed to their solar value. Other parameters 
were all free; the electron temperature kT c , the ionization 
timescale r (a product of the electron density and the 
elapsed time after the shock heating), and the emission 
measure (= / n e n*\dl, where n e and tir are the number 
densities of hydrogen and electrons and dl is the plasma 
depth) . We also allowed the column density Nn to vemis- 
sion measureary freely. For the hard X-ray tail Katsuda 
et al. (2009) applied the CXB model plus an additional 
power law component. They used the same value of N h 
for the CXB as for other components. In our analysis, it 
is fixed to 3.6xl0 21 cm -2 calculated by nh FTOOL 1 and 
separated from the Nu for G156.2+5.7. The additional 
power law component was not applied for our analysis 
unless the CXB surface brightness exceeded the literature 
value. 

Katsuda et al. (2009) found excess emission over the 
CXB above 3 keV in the NW-rim and argued that the ob- 
served excess originates from a non- thermal emission from 
G156.2+5.7. In their analysis, they excluded a region of 
1'.2 in radius where 2XMM J045637.2+522411 is located. 
Although Yamauchi et al. (2010) estimated the exten- 
sion of the X-ray emission from 2XMM J045637.2+522411 
to be at least ~50" (the surface brightness within it is 
~2.64x 10“ 10 erg cm~ 2 s -1 deg~ 2 in the 2.0-10.0 keV band) 
based on the XMM-Newton data, we noticed that its ex- 
tent might be larger. Figure 3 left and right panels show 
a full-band (0.3-8. 0 keV) and a hard-band (2. 0-8.0 kev) 
image, respectively. In figure 3, we show the small cir- 
cular region of 1'.2 in radius as used by Katsuda et al. 
(2009). We find that the X-ray emission from the clus- 
ter extends beyond this circle even in the hard band. 
Figure 4 shows a radial profile from the center of 2XMM 
J 045637.2+5224 11 in the 2. 0-8.0 keY band. The dotted 

http://heasarc.gsfc.nasa.gov/cgi-bin/Tools/w3nh/w3nh.pl 



Fig. 3. Full-band (0.3-8.0 keV; left), and hard-band (2.0-8.0 
keV; right) images of NW_rim. The FbV of the XIS axe shown 
with white squares. The radii of smril and large circles are 
1/.2 and 4 , .0, respectively. 
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Fig. 4. Radial profile of NW_rim in the 2.0-8.0 keV energy 
band. The center corresponds to the position of 2XMM 
J 045637. 2+522411. The dotted vertical line represents the 
radius of 1'.2. 

vertical line corresponds to the radius of l 7 . 2. Figure 4 
shows that the surface brightness decreases to half at the 
radius of ~3'. To rule out the possible contamination from 
2XMM J045637.2+522411, we therefore excluded a circu- 
lar region with radius 4'.0 as shown with the white large 
circle in figure 3. The fitting result after excluding this 
larger region is shown in figure 5. The best-fit parameters 
are shown in table 3. 

The surface brightness of the hard component in the 
NW_rim is calculated to be 1.88xl0~ n erg cm” 2 s“ 1 deg~ 2 
(2.0-10.0 keV; table 3). On the other hand, the CXB 
surface brightness is estimated to be 1.57±0.32x 10^ n erg 
cm“ 2 s _1 deg~ 2 based on blank sky observations with the 
XIS (Tawa 2008). In their analysis, a spatial fluctuation 
14% (1 <t) is considered for the XIS FoV (Tawa et al. 2008). 
Accordingly, the surface brightness of the hard component 
in the NW-rim (and also E_rim and S-rim) is within the 
uncertainty of the CXB surface brightness. We conclude 
that no additional power law component is required for 
our fit. 

3.3. Model Fit for the Inner Regions' Spectra 

We also applied the single-component VNEI model for 
each inner region. However, this model was inadequate 
to fit the data due to its simplicity. As previously shown 
by Katsuda et al. (2009), the thermal emission from the 
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No. 


Component 

Parameter 


Value 




NW_rim 

E_rim 

S_rim 

CXB model 





absorption 

Nh [10 21 cm -2 ] . . . 


3.6 (fixed) 


Broken Power-law 

Ti 


2.0 (fixed) 



Ebreak [keV] 


0.7 (fixed) 



r 2 


1.4 (fixed) 



surface brightness* 

1.88xl0~ u 

1.61 xlO -11 

1.56X10" 11 

LHB model 





APEC 

kT e [keV] 


0.1 (fixed) 



Abundances 


1 (fixed) 



surface brightness* 

3.3 xlO -16 (fixed) 


VNEI model 





absorption 

Nh [10 21 cm -2 ] . . . 

2.3+0.1 

2. 0+0.1 

3. 1+0.1 


kT e [keV] 

0.45+0.01 

0.49+0.01 

0.41+0.01 


log(r [s cm -3 ]) . . . 

10.64+0.02 

10. 54^02 

10.80l°;°3 


N 

0.12+0.01 

0.15+0.01 

0.50±0.01 


0 

0.24+0.01 

0.22+0.01 

0.25 ±0.01 


Ne 

0.42+0.01 

0.44+0.01 

0.57±0.01 


Mg 

0.37+0.02 

0.34+0.01 

0.47±0.01 


Si (=S) 

0.71+0.08 

0.32+0.01 

0.48±0.14 


Fe (=Ni) 

0.36+0.01 

0.37+0.01 

0.40±0.01 

^ O n 1 n lr/\A 7 rv-w-n -wfrx , 

X 2 /dof 

~ Ui. ■ n. r 

1395/1245 

1441/1190 

T 

1183/1031 


Other elements are fixed to solar values. 

The errors are in the range Ax 2 < 2.7 on one parameter. 
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Fig. 5. NW_rim spectrum fitted with the single-component 
VNEI model. 

interior of G 156.2+5. 7 consists of two components with 
different origins; a high-temperature ejecta component 
and a low-temperature ISM component. Thus we applied 
a two-component VNEI model for the spectra obtained 
from Center, East and NW. In this model, we fixed the 
ISM abundances to the average values based on the fit- 
ting results of rim observations (see table 3). The av- 
erage abundances of the ISM component are set as fol- 
lows: N=0.23, 0=0.24, Ne=0.45, Mg=0.39, Si(=S)=0.50, 
Fe(=Ni)=0.38. We fixed the other abundances to the solar 
values. We also set kT e to be 0.45 keV which is an average 
electron temperature for the rim regions of G156.2+5.7. 
The values of r and emission measure were allowed to vary. 


Meanwhile, in the ejecta component, the abundances of 
O, Ne, Mg, Si, S and Fe were free while we set the abun- 
dances of C and N equal to O, Ni equal to Fe. All other 
abundances were fixed to their solar values. The other 
free parameters were kT ej r, emission measure, and the 
column density Nft. 

We note that the temperature of the low-kT component 
may be higher than that of the rim because the line-of- 
sight ISM component includes ambient medium shocked 
earlier in the evolution of the remnant when the shock 
velocity was higher. In general, higher ISM temperatures 
will increase the line emission from species such as Si and 
S, which may affect the ejecta abundance measurement. 
Since it is impossible to measure the real average tem- 
perature of the ISM component, we examined how the 
metal abundances change by increasing the ISM temper- 
ature from 0.45 keV to 2.0 keV. As a result, we concluded 
that the ISM temperature hardly affect the ejecta abun- 
dances within the range of the error. We thus fixed the 
ISM temperature to be 0.45 keV for the inner regions’ 
spectra. 

4. Discussion 

4-1. Ejecta Structure 

Figure 6 shows band images of Ovil, Ovill, NelX, 
Mg XI, Si XIII and S XV using all the Suzaku data. In figure 
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Fig. 6. Band images of several elements. The white contours were derived from a soft X-ray image taken by ROSAT PSPC. 


6. both the Si and the S distributions appear to be cen- 
trally concentrated, while those of 0 VII and 0 VIII show 
a clear limb brightening, which suggests that the observed 
heavy elements originated from the ejecta. Furthermore, 
figure 6 suggests that the distribution centers of Si and S 
shift toward the northwest. 

To examine the radial distribution of the plasma, we 
divided the data into annular regions, as shown in fig- 
ure 1 right. As the annular center we used 04 h 58 m 54 s .9, 
51°44 / 36".l (J2000), the optical axis during the observa- 
tion of Center. We used an annular width of 2', consis- 
tent with the angular resolution of the XIS. In order to 
estimate the boundary of the ejecta-dominated region, we 
applied both the single-component VNEI model and the 
two-component VNEI model to the spectra obtained from 
the E_rim and NWjrim. We determined which model is 
preferred by using the F-test with a significance level of 
99%. We found that the second component is required 
in the NW_rim only within roughly ~85% of the shock 
radius (~52'). Figure 7 shows example spectra obtained 
from the annular regions 25' away from the annular cen- 
ter toward the east (E-25; left) and the northwest (NW- 
25; right). The best-fit parameters are shown in table 
4. While both spectra require the two-component VNEI 
model, the abundances of the ejecta component are signifi- 
cantly different particularly for Si and S. We also note that 
there is no positive evidence that G156.2+5.7 is an O-rich 
SNR as claimed by Hudaverdi et al. (2010). Their argu- 
ment is that the south and east side of their FoV are soft 
X-ray dominated (0. 4-1.0 keV), suggesting the abundant 
presence of 0. However, we consider that such excesses 


are mainly due to the contribution of the low-temperature 
ISM component. 
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Fig. 7. Example spectra of two annular regions in 
G156.2+5.7. Both were fitted with the two-component VNEI 
model. 

Figure 8 shows the radial profiles of the emission mea- 
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Table 4. Spectral parameters (see figure 7) 
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Component 

Parameter 

Value 



Region E-25 Region NW-25 

CXB model 



absorption 

Nr [10 2i cm -2 ] 

3.6 (fixed) 

Broken Power-law 

r x 

2.0 (fixed) 


Fbreak [keV] 

0.7 (fixed) 


r 2 

1.4 (fixed) 


surface brightness* 

1.48xl0 -11 1.74xl0~ n 

LHB model 



APEC 

kT e [keV] 

0.1 (fixed) 


Abundances 

1 (fixed) 


surface brightness* 

3.3xl0- 16 (fixed) 

VNEI model 



absorption 

Nh [10 21 cm -2 ] 

1.8±0.1 3.4±0.1 

Low-kT e component 

kT e [keV] 

0.45 (fixed) 


log(r [s cm -3 ]) 

11.381“;°? 10.001°;°° 


N 

0.23 (fixed) 


0 

0.24 (fixed) 


Ne 

0.45 (fixed) 


Mg 

0.39 (fixed) 


Si 

0.50 (fixed) 


S 

0.50 (fixed) 


Fe (=Ni) 

0.38 (fixed) 

High-kT e component 

kT e [keV] 

0-59lgS? 0.58±0.01 


log(r [s cm -3 ]) 

11.29±0.04 11.14lgoJ 


O (=C=N) . 

118lo2i 0.61±0.09 


Ne 

0.51±0.11 0.56±0.08 


Mg 

0.83lg 0.63±0.08 


Si 

0.48±0.22 2.25±0.20 


S 

<1.5 3.40±0.96 


Fe (=Ni) 

1.05±0.07 1.45±0.06 

* on inn 1 1 

X 2 / dof 

379/326 533/453 


Other elements are fixed to solar values. 

The errors are in the range Ax 2 < 2.7 on one parameter. 


sure for various metals (O, Ne, Mg, Si, S, and Fe) in the 
ejecta. We found prominent central concentrations for Si, 
S, and Fe, while the lighter elements (O, Ne, and Mg) are 
distributed more uniformly. The emission measure of Si 
peaks ^lO 7 from the annular center and the overall dis- 
tribution is quite similar to that of S. The co-existence 
between Si and S is consistent with a theoretical picture 
that they are mainly produced in the static O-burning 
layer during the stellar evolution and the incomplete Si- 
burning layer during the SN explosion. Similar shifts are 
observed in other evolved SNR s, for example, the Cygnus 
Loop (Uchida et al. 2009). If such inhomogeneities are 
real, one of the possible causes is an asymmetric explosion 
of a progenitor star. While many SN models have been 
proposed in recent decades, a model invoking delayed ex- 
plosion caused by neutrino heating is worth considering in 
the context of these data. For example, Marek &: Janka 
(2009) simulated an SN explosion of 15M© caused by the 


SASI-aided core-collapse (SASI; standing accretion shock 
instability). They expect that the heavy elements, such 
as Si and Fe, are ejected in one direction, while the lighter 
elements, such as O, Ne, and Mg, are ejected in all di- 
rections. However, since our FoV is limited, we cannot 
determine the precise distribution center of each ejected 
element. Also, the annular center does not necessarily 
correspond to the explosion center. It should be care- 
fully studied whether the origin of such inhomogeneous 
distribution is actually the result of an asymmetric SN 
explosion or not. 

4.2. Origin of hard X-ray emission 

As shown in section 3.2, we found no robust evidence 
for non-thermal emission in the NW_rim when we ex- 
cluded the region contaminated by the cluster 2XMM 
J045637.2+522411. The consistency of the measured 
CXB surface brightness from the NW-rim with other ob- 
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Fig. 8. Distributions of emission measure for various metals (O, Ne, Mg, Si, S, and Fe) in the ejecta. Open circle, cross, square, 
and triangle correspond to the spectral extracted regions, East, Center, NW, and NW_rim, respectively. The horizontal axis shows 
a distance from the center of Center. The positive and the negative values represent directions toward the northwest and the east, 
respectively. 
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servations’ measurements (see table 2) as well as with val- 
ues found in the literature strengthens the conclusion that 
there is no non-thermal emission from the NW_rim. 

We next investigate spatial fluctuations of the hard X- 
ray surface brightness by using all the spectra obtained 
from the annular regions described above (see figure 2). 
Prom figure 2 left, there is no correlation between the 
radio synchrotron and the hard X-ray emission contrary 
to the prediction by Katsuda et al. (2009). The hard 
X-ray distribution does not have a bilaterally symmetric 
morphology like some non-thermal SNRs (e.g., SN1006; 
Koyama et al. 1995). Figure 2 right shows a histogram of 
the calculated surface brightnesses. The horizontal arrow 
represents the 1 -a spatial fluctuation of the CXB surface 
brightness estimated by Tawa (2008). Surface bright- 
nesses in almost all regions (~94%) are within the 3-cr 
fluctuation, which suggests that the observed hard X-ray 
emission is dominated by the CXB, Although Katsuda 
et al. (2009) argued that non-thermal emission was also 
detected at the SNR’s center, figure 2 shows that the hard 
surface brightness from this region is also consistent with 
CXB fluctuations, as is the case in the NW«rim. As a re- 
sult, there is no robust evidence for non-thermal emission 
from this remnant. 


4.3. Estimating the age of G156.2+5.7 


As shown above, the observed hard X-ray emission is 
dominated by the CXB everywhere except in the region of 
4'.0 around 2XMM J045637.2+522411. This fact supports 
the idea that G156.2+5.7 is an evolved SNR, despite the 
claim of Gerardy & Fesen (2007) that it was a young 
nearby one. As explained in section 4.1, G156.2+5.7 has 
a layered metallicity structure that is observed in other 
evolved SNRs (e.g., Cygnus Loop; Tsunemi et al. 2007). 
Such structure suggests that the reverse shock has already 
reached at the center of G156.2-f-5.7. Furthermore, the 
electron temperature (< 1 keV) is much lower than that 
in many young SNRs. All the results support the picture 
of an evolved SNR. For these reasons, we prefer ~1.1 kpc 
(Pfeffermann et al. 1991) to ~0.3 kpc (Gerardy & Fesen 
2007) for its distance d in the following analysis. 

Assuming thermal equilibrium between electrons and 
protons, we find that the forward shock velocity v 3 is esti- 
mated to be v 3 ~ 614(fcT e /0.45keV) 0,5 km s -1 . However, 
Ghavamian et al. (2001) calculate an initial equilibration 
of only 50% for shock speeds of ~ 600km s -1 . In such a 
non-equilbrium case, the electron temperature kT e , and 
the proton temperature, kT p are expressed in the follow- 
ing equations (Ghavamian et al. 2001): 


K-./eq + (2 - /eq)m p ] v ] , 

kTe ~\h> [ m ^ eq + ( 2 “ u s> 
Ik — 

Tp~ 2 -/ eq ’ 


( 1 ) 

( 2 ) 

( 3 ) 


where m p>e are the proton and electron masses, respec- 
tively. If the initial electron-proton equilibration, / eq is 
0.5, the forward shock velocity is 960(fcT e /0.45keV)° 5 


km s -1 , which gives an upper limit of v s . 

On the other hand, the ambient densi- 

ties (no) for the NW-rim and the E-rim are 
no,NW = 0.084(d/l.l kpc) -1 cm -3 and no,E = 

0.091 (d/1.1 kpc) -1 cm” 3 , respectively, by fitting their 
emission measure profiles with a Sedov model (Katsuda 
et al. 2009). Due to the lack of statistics, it is hard to 
estimate the ambient density for the S_rim in the same 
way. We therefore consider the average value of two cal- 
culated no to be a global ambient density of G156.2-I-5.7: 
hence n 0 = 0.088(d/l.l kpc) -1 cm -3 . Applying a simple 
Sedov analysis, the age £4 = i/10 4 vr is as follows: 


088 cm -3 


)■ 


„( V, N-5/3 / Eo y'*( 

4 V614kms-V \10 51 erg/ W 
for the thermal equilibrium case, and 

, .. 0 -/ y*t*( Eo \ 1/3 ( no \ 

4 VOfiOkms -1 / \10 51 ergy \0.088cm -3 / 


for the non-equilibrium case. Here E 0 is the explosion 
energy. As a result, we prefer the estimate from the 
ASCA observation (15000 yr; Yamauchi et al. 1999) to 
that from the ROSAT observation (26000 yr; Pfeffermann 
et al. 1991). In any cases, G156.2+5.7 is not a young SNR, 
as claimed by Gerardy & Fesen (2007). While Katsuda 
et al. (2009) pointed out that non-thermal X-ray emission 
would be detectable even from such an evolved SNR if 
the flux ratio of non-thermal to thermal emission is higher 
than those in a typical evolved SNR, it is more reasonable 
to conclude that G 156. 2+5. 7 is, itself, a typical evolved 
SNR, since CXB exceeds that of any X-ray synchrotron 
component. 


5. Summary 

We performed a set of 6 pointing observations of 
G 156.2+5.7 with Suzaku. The spectra were well fit- 
ted with two-component or single-component VNEI (with 
CXB and LHB) models for the inner regions and rim re- 
gions, respectively. Our results are consistent with the 
conclusion of Katsuda et al. (2009), that these compo- 
nents represent reverse shock-heated ejecta and forward 
shock-heated ISM, respectively. 

From the ejecta component, we found prominent central 
concentrations of Si, S, and Fe, while the lighter elements 
(O, Ne, and Mg) are distributed more uniformly. Such 
distributions reflect a layered-metallicity structure of the 
progenitor star. 

A single-component VNEI model provided an accept- 
able fit for all the rim regions; no additional power law 
component was required for our fit. We found no ro- 
bust evidence for the non-thermal emission in the NW-rim 
when we excluded the region contaminated by the cluster 
of galaxies, 2XMM J045637. 2+522411. Accordingly, we 
conclude that the hard X-ray emission in G156.2+5.7 is 
sufficiently explained as CXB emission. 

The estimated forward shock velocity is relatively 
slow compared with young SNRs: v 8 ~ 614 — 
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960 (kT e /Q 45keV) 0 5 km s” 1 . We estimated its age to Yamauchi, S., Koyama, K., Tomida H., Yokogawa J., & 
be 7,000-15,000(d/l.lkpc) _1 yr; therefore, G156.2+5.7 is Tamura K. 1999, PASJ, 51, 13 

not a young, but an evolved SNR. Yamauchi, S., Ueno, M., Bamba, A. s & Koyama, K. 2010, 

PASJ, 62, 219 

Acknowledgments Wilms, J., Allen, A., & McCray, R. 2000, ApJ, 542, 914 

H.U. is supported by Japan Society for the Promotion of 
Science (JSPS) Research Fellowship for Young Scientists. 

H.Y. is supported by JSPS Research Fellowship for 
Research Abroad. The work is partially supported by 
the Ministry of Education, Culture, Sports, Science and 
Technology (Grant-in- Aid No.23000004). 

References 

Anders, E., & Grevesse, N. 1989, Geochim. Cosmochim. Acta, 

53, 197 

Amaud, K. A. 1996, in ASP Conf. Ser. 101, Astronomical 
Data Analysis Software and Systems V, ed. G. H. Jacoby 
& J. Barnes (San Francisco: ASP), 17 
Borkowski, K. J., Lyerly, W. J., & Reynolds, S. P. 2001, ApJ, 

548, 820 

Enoto, T., et al. 2009, ApJL, 693, L122 

Gerardy, C L., & Fesen R. A. 2007, hlNRAS, 376, 929 

Ghavamian P., Raymond, J., Smith, R. C., & Hartigan, P. 

2001, ApJ, 547, 995 

Hudaverdi, M., Bozkurt, M., Sezer, A., Gok, F., & Ercan, 

E. N. 2010, 9th International Conference of the Hellenic 
Astronomical Society, 424, 179 
Ishisaki, Y., et al. 2007, PASJ, 59, SI 13 
Kaplan, D. L., Gaensler. B. M., Kulkarni, S. R., & Slane, P. O. 

2006, ApJS, 163, 344 

Katsuda, S , Petre, R., Hwang, U., Yamaguchi, H., Mori, K., 

& Tsunemi, H., 2009, PASJ, 61, S155 
Kawaguchi, T., Nakanishi, K. Kohno, K., Ohta, K., & Aoki, 

K., 2008, ApJ, 676, 137 

Kesteven M. J., & Caswell, J. L., 1987, A&A, 183, 118 
Koyama, K, Petre, R., Gotthelf, E. V., Hwang, U., Matsura, 

M., Ozaki, M., & Holt, S. S., 1995, Nature, 378, 255 
Koyama, K et al. 2007, PASJ, 59, S23 
Lazendic, J. S., & Slane, P. O. 2006, ApJ, 647, 350 
Marek, A. & Janka, H.-Th. 2009, ApJ, 694, 664 
Miller, E. D., et al. 2008, PASJ, 60, S95 
Pfeffermann E., Aschenbach B., & Predehl P. 1991, A&A, 

246, L28 

Reich W., Fuerst E., k Amal E. M. 1992, A&A, 256, 214 
Revnivtsev. M., Sazonov, S., Gilfanov, M., Churazov, E., k 
Sunyaev, R. 2006, A&A, 452, 169 
Smith, R. K., Brickhouse, N. S., Liedahl, D. A., & Raymond, 

J. C. 2001, ApJL, 556, L91 

Sugizaki, M., Mitsuda, K., Kaneda, H., Matsuzaki, K., 

Yamaguchi, S., & Koyama K. 2001, ApJ, 134, 77 
Tawa, N., et al. 2008, PASJ, 60, Sll 
Tawa, N. 2008, PhD thesis, Osaka University 
Tsunemi, H., Katsuda, S., Norbert, N., & Miller, E. D. 2007, 

ApJ, 671, 1717 

Uchida, H., Tsunemi, H., Katsuda, S. Kimura, M., & Kosugi, 

H. 2009, PASJ, 61, 301 

Xu, J. W., Han, J. L., Sun, X. H., Reich, W., Xiao, L., Reich, 

P., Wielebinski, R. 2007, A&A, 470, 969 
Yamauchi, 3., Ueno, S., Koyama, K., Nomoto, S., Hayashida, 

K. , Tsunemi, H., Asaoka, I. 1993, PASJ, 45, 795 



